Gamma rays and neutrinos from the Crab Nebula produced by pulsar accelerated nuclei 
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We investigate the consequences of the acceleration of heavy nuclei (e.g. iron nuclei) by the 
Crab pulsar. Accelerated nuclei can photodisintegrate in collisions with soft photons produced in 
the pulsar's outer gap, injecting energetic neutrons which decay either inside or outside the Crab 
Nebula. The protons from neutron decay inside the nebula are trapped by the Crab Nebula magnetic 
field, and accumulate inside the nebula producing gamma-rays and neutrinos in collisions with the 
matter in the nebula. Neutrons decaying outside the Crab Nebula contribute to the Galactic cosmic 
rays. We compute the expected fluxes of gamma-rays and neutrinos, and find that our model could 
account for the observed emission at high energies and may be tested by searching for high energy 
5— I , neutrinos with future neutrino telescopes currently in the design stage. 

PACS numbers: 98.70.Rz, 98.58.Mj, 95.85.Ry, 97.60.Gb, 98.70.Sa. 
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The Crab Nebula is a well established 7-ray source 
~H ■ with a complex spectrum extending up to at least TeV 
energies |lj . Its emission below a few tens of MeV is inter- 
preted as synchrotron emission by electrons with energies 
up to ~ 10 15 ~ 16 eV in the Crab Nebula magnetic field of 
strength ~3x 10~ 4 G ||. These models interpret the 
emission at higher energies as being due this same papu- 
la , lation of electrons by inverse Compton scattering (ICS) 
O'S ' of synchrotron photons (SSC models). The TeV 7-ray 
emission might also originate near (but above) the Crab 
O { pulsar light cylinder as a result of ICS scattering of in- 
frared photons produced in one of the pulsar's outer gaps 
by pairs escaping from another outer gap ||. The 
possible importance of hadronic processes as sources of 
7-rays inside the Crab Nebula has also been considered 
by Cheng et al. who calculated the high energy 7-ray 
. £h ! fluxes assuming that relativistic protons accelerated in 
the Crab pulsar outer gap interact with the matter inside 
the Nebula, and note that this process may contribute in 
. . . 1 the TeV 7-ray range. 

In this Letter, we analyze the consequences of acceler- 
ation of heavy nuclei in the pulsar magnetosphere as a 
possible mechanism of energetic radiation from the Crab 
Nebula. The importance of photodisintegration of nuclei 
is determined by the reciprocal mean free path which we 
calculate using cross sections of Karakula & Tkaczyk || 
and show in Fig. [l] for nuclei with atomic mass numbers 
between 4 and 56 during propagation in the radiation 
field supposed to be present in the outer gap of the Crab 
pulsar. Ho Q has shown that the total luminosity and 
spectrum of the outer gap radiation is not strongly de- 
pendent on P. For P = 33 ms we adopt the photon 
number density given by Eq. 7.1 in Paper II of Cheng et 
al. 0, and for other periods we scale this by r^c- We 
assume this radiation is produced by cascading due to 
electrons and positrons accelerated in the gap, and so is 
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highly anisotropic and directed across the gap (in both 
directions). The dimension of the outer gap is of the or- 
der of the radius of the light cylinder, ^ gap w tlc = c/fl. 
For the case of the Crab pulsar, tlc ~ 1-5 X 10 8 cm, and 
it is evident from Fig. [j] that multiple photodisintegra- 
tions of primary Fe nuclei will occur provided they can 
be accelerated to Lorentz factors above ~ 10 3 . To check 
whether pion photoproduction and e ± pair production 
by Fe nuclei are important energy loss mechanisms in the 
outer gap, we compute also the mean interaction length 
for pion photoproduction using the cross section given 
by Stecker Q, and the energy loss length for pair pro- 
duction using the approximations given by Chodorowski 
et al. ||. The results are shown in Fig. [|. It is clear 
that these processes are not important loss mechanisms 
compared to photodisintegration, e.g the pair production 
energy loss length is two orders of magnitude longer than 
the length of the gap. 

We assume that Fe nuclei can escape from the polar 
cap surface of the Crab pulsar, and move along magnetic 
field lines to enter the outer gap where they can be ac- 
celerated in the outer gap potential as in the model of 
Cheng et al. and Ho ||. The nuclei accelerated in 
the outer gap will interact with photons either produc- 
ing secondary e pairs (with negligible loss of energy) or 
extracting a nucleon. The pairs could eventually saturate 
the electric field of the gap although this is far from cer- 
tain. We shall assume that for some reason the field is not 
shorted. Possibilities include: (a) nuclei, having much 
larger Larmor radii and much lower synchrotron and cur- 
vature losses than electrons, more easily drift across the 
magnetic field lines away from the region of the gap where 
they had produced pairs; (b) the secondary pairs may 
not be able to saturate the electric field in the outer gap 
which is close to the light cylinder since the Goldreich 
& Julian jlO| density tends to infinity where the gap ap- 
proaches the light cylinder. 

In order to obtain the energy spectrum of neutrons 
extracted from Fe nuclei, iV n (7 n ), we simulate their ac- 
celeration and propagation through the outer gap using 
a Monte Carlo method. To obtain the rate of injection 
of neutrons per unit energy we multiply N n {^ n ) by the 
number of Fe nuclei injected per unit time, Np e , which 
can be simply related to the total power output of the 
pulsar L Crah (B,P) [0, 

NFc = £L Crah (B, P)/Z*(B, P), (1) 

where £ is the parameter describing the fraction of the 
total power taken by relativistic nuclei accelerated in the 
outer gap, Z = 26 is the atomic number of Fe, B is the 
surface magnetic field, P is the pulsar's period, and 

is the potential difference across the outer gap We 
assume B = 4 x 10 12 G. 
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Soon after the supernova explosion, when the nebula 
was relatively small, nearly all energetic neutrons would 
be expected to decay outside the nebula. However, at 
early times we must take account of collisions with mat- 
ter. The optical depth may be estimated from r n H ~ 
a pp n H r w 8.6 x 10 14 MiV$ i~ 2 , where M = M\Mq is 
the mass ejected during the Crab supernova explosion in 
units of solar masses, r = vt, uh — M/(4/3irr 3 m p ) is 
the number density of target nuclei, and v = 10 8 vg cm 
s _1 is the expansion velocity of the nebula. We note that 
r nH = 1 at t » 0.93M 1 1/2 w^" 1 y. 

The number and spectrum of relativistic protons from 
neutron decay at the present time is determinated by 
the injection rate of neutrons into the Crab Nebula in- 
tegrated over time since the pulsar's birth. We estimate 
the evolution of the Crab pulsar's period from birth to 
the present time taking account of magnetic dipole radi- 
ation energy losses and gravitational energy losses for an 
ellipticity of 3 x 10~ 4 Magnetic dipole losses deter- 
mine the pulsar period at present, but the initial period 
is determined largely by gravitational losses and is proba- 
bly shorter than ~ 10 ms. Hence, we consider two initial 
periods: 5 ms, and 10 ms. 

The spectrum of protons from neutron decay outside 
the Crab Nebula is given by 

N° ut (l P ,t C N) = 

""* dtN Fc (t)N n ( lp ,t)e- T - H ^e- vt ^/ c ^ (3) 

where r„ rj 900 seconds is the neutron decay time, and 
we make the approximation that the Lorentz factor of 
protons is equal to that of parent neutrons, 7 P w j n . 

The equation for the spectrum of protons injected in- 
side the Crab Nebula is more complicated because we 
must take account of proton-proton collisions and adia- 
batic energy losses due to the expansion of the nebula. 
The Lorentz factor of these protons at time t after the 
explosion such that their present Lorentz factor is 7 P is 
given by 

where T pp (t) is the optical depth for collision of pro- 
tons with matter between t and £cn, and is given by 
T pp (t) w 1.3 x 10 17 Mi«^ 3 (t- 2 -^), and K rj 0.5 is the 
inelasticity coefficient in proton-proton collisions. 

Since we are interested in protons interacting inside the 
nebula at the present time, we must also include those 
neutrons which decayed at locations outside the nebula 
at time t but which will be inside the nebula at time 
icN- Taking account of all these effects, we arrive at the 
formula below for the proton spectrum inside the nebula 
at time ten, 

/■*CN 

1 / MN„J+\f>- T nH{t) 



N p n ( lpi tax) = % / dtN Fe (t)e 
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x N n (^ p (t),t)j p (t)[l - exp(-vt/cry p (t)r n )] 



+ 




dt'N n {j p (t'),t) lp (t') 



v exp(-vt' / cj p (t')T n ) 

C7p(*') r n 



(5) 



The first term gives the contribution from neutrons de- 
caying initially inside the nebula while second term gives 
the contribution from neutrons decaying at points ini- 
tially outside the nebula which will be inside the nebula 
at time icN- 

We compute separately the spectra of protons from 
neutrons decaying inside and outside the Crab Nebula. 
In Fig. ^ proton spectra are shown for two initial peri- 
ods and present nebula radii of 1 pc and 2 pc. To check 
if protons decaying inside will remain near where they 
were produced, or diffuse into the Galactic environment, 
we estimate their typical diffusion distance during time 
*CN, Xdif ~ (crf^cN/3) 1 / 2 , where tl = m p j p c 2 /eB is the 
Larmor radius, and we have used the minimum diffusion 
coefficient. For protons at the peak of the energy spec- 
trum with 7 p = 10 5 in a magnetic field of the order of 
B = 5 x 10~ 6 G, the typical diffusion distance is Xdif ~ 1 
pc, which is comparable to the radius of the Crab Neb- 
ula. Hence protons which decayed inside, were proba- 
bly trapped and accumulated inside the nebula. Protons 
from neutrons decaying outside the nebula will typically 
have higher energies and diffuse farther, thus escaping to 
become Galactic cosmic rays. 

The protons which have accumulated inside the Crab 
Nebula since the pulsar was formed can produce observ- 
able fluxes of 7-rays and neutrinos. We compute the 
expected 7-ray spectra for four different models, taking 
various initial pulsar periods, present sizes and masses of 
the Crab Nebula. The fluxes may possibly be enhanced if 
protons are efficiently trapped by the dense filaments as 
suggested by Atoyan and Aharonian [fij| . Filaments with 
density ~ 500 cm -3 are present in the Crab Nebula [fTi) . 
Therefore we introduce an effective density experienced 
by the protons inside nebula, n c ^ = fina, where tih is 
defined above, and the parameter /i takes into account 
the possible effects of proton trapping by the filaments. 
The 7-ray spectra are computed in terms of a scaling 
model , and the photon fluxes expected at the Earth 
are shown in Fig. ^| for £/i = 1, and are compared with 
observations of the Crab Nebula above 0.2 TeV. Results 
are shown for models I to IV having Pq, tcn and Mi 
as specified in Table EL and for a distance to the Crab 
Nebula of 1830 pc flj. Comparison with the Whipple 
observations at 10 TeV allows us to place constraints on 
the free parameters of the model, and upper limits on £/i 
which are given in Table Q as (£/i) 7 . Note that it is usu- 
ally argued for the standard pulsar model that the rate 
of injection of Fe nuclei into the pulsar magnetosphere 
should not cause the charge density to exceed the Gol- 
dreich & Julian density |l(J . In the case of no additional 
currents flowing through the magnetosphere, this condi- 
tion constrains the value of £ to £ « 0.1 for the Crab 
pulsar with period 33 ms. However, in the presence of 
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additional currents the value of £ can be higher. 

The question of the importance of hadronic interac- 
tions in the Crab Nebula can be settled by the detection 
of a neutrino signal from the Crab. In Fig. ^ we show the 
neutrino spectrum produced in collisions of protons with 
matter inside the Crab nebula for the models considered 
above, and compare this with the atmospheric neutrino 
background flux within 1° of the source direction ^0|. It 
is clear that neutrino detectors with good angular reso- 
lution should be able to detect neutrinos at 10 TeV from 
the Crab Nebula if £/i is greater than given in Ta- 

ble H Note that in all cases (£fJ,) v < (£^) 7 and so the 
possibility of v detection is allowed by the existing 7-ray 
observations. 

In the present paper we have not considered those nu- 
clei which had survived propagation through the outer 
gap radiation field and were injected into the inner neb- 
ula region following magnetic field lines. These nuclei 
(and protons) are expected to accumulate in the inner 
part of the Crab Nebula were the magnetic fields are 
^2xlO~ 3 G measured in the synchrotron wisps which 
are located at ~ 0.2 pc The diffusion distance scale 
during the age of the Crab Nebula is ~ 0.3 pc for a 
Lorentz factor 7^4 = 10 6 and B — 10~ 3 G. This means 
that most nuclei will probably be confined to the inner 
part of the Crab Nebula where there is no evidence of 
dense matter H] (other than the neutron star). Hence, 
we may be justified in neglecting the contribution of such 
nuclei to the 7-ray and neutrino production by hadronic 
collisions in comparison with the contribution from in- 
teractions of protons from neutron decay. However, some 
fraction of the energy of these nuclei could be transferred 
to very high energy electrons by resonant scattering |22j . 
Such electrons might then produce additional 7-rays by 
the SSC process as in the models discussed in ref. 
Hence, we suggest that the complex Crab Nebula spec- 
trum may in fact be formed as a result of both processes, 
i.e. SSC and as a result of pulsar acceleration of nuclei 
discussed here. 

Another consequence of our model is that particles (nu- 
clei and protons from neutron decay) will be also injected 
by the Crab Pulsar and other pulsars into the Galac- 
tic cosmic rays, this may effect the cosmic ray spectrum 
and composition. This problem needs closer investiga- 
tion, and should take into account the parameters of the 
Galactic population of pulsars at birth. However, it is 
worth noting that the energy distribution of the protons 
from decay of neutrons outside the Crab nebula peaks at 
~ 10 15 eV (see Fig. ||), which is close to the energy of the 
knee in the cosmic ray spectrum. 
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(£/i) 7 at 10 TeV 
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FIG. 1. Reciprocal mean free paths for photodisintegra- 
tion of nuclei with mass numbers, A = 56, 40, 32, 24, 16, 8, 
and 4 (full curves from the top) in the radiation field of the 
Crab pulsar's outer gap as a function of the Lorentz factor of 
the nuclei. The reciprocal mean free path for pion photopro- 
duction (dotted curve) and the energy loss length for e ± pair 
production (dashed curve) by iron nuclei are also shown. 
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FIG. 2. Spectra of protons from the decay of neutrons in- 
jected by the Crab pulsar inside and outside the Crab Nebula 
assuming Mi = 3 for the following parameters: Po = 10 ms, 
and Crab Nebula radius rcN = 1 pc (dot-dashed curve) and 
2 pc (dashed curve), and for P — 5 ms and rcN = 1 pc (full 
curve) . 
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FIG. 3. Spectra of 7-rays from interactions of protons with 
matter inside the Crab Nebula for the different proton spectra 
shown in Fig ^| and for the different models of the nebula (I 
- IV) considered in the text. Observations: Whipple Obser- 
vatory Jl6| (dotted line and error box); THEMISTOCLE @ 
(+); and CANGAROO (is| (solid line and error boxVUpper 
limits from various experiments mentioned in ref. |19(: T - 
Tibet, H - HEGRA, C - CYGNUS, and U - CASA-MIA. 
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FIG. 4. Spectra of neutrinos (y^ + Vfj) expected for the dif- 
ferent models (I - IV) considered in the text. The atmospheric 
neutrino background within 1° of the source is indicated 
by the hatched band which shows the variation with angle: 
horizontal (upper bound), vertical (lower bound). 
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